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In this paper, the gravitational field equations for static spherically sym- 
metric perfect fluid models with a polytropic equation of state, p = fep^+'^Z", 
are recast into two complementary 3-dimensional regular systems of ordinary 
differential equations on compact state spaces. The systems are analyzed nu- 
merically and qualitatively, using the theory of dynamical systems. Certain 
key solutions are shown to form building blocks which, to a large extent, de- 
termine the remaining solution structure. In one formulation, there exists a 
monotone function that forces the general relativistic solutions towards a part 
of the boundary of the state space that corresponds to the low pressure limit. 
The solutions on this boundary describe Newtonian models and thus the rela- 
tionship to the Newtonian solution space is clearly displayed. It is numerically 
demonstrated that general relativistic models have finite radii when the poly- 
tropic index n satisfies < n ^ 3.339 and infinite radii when n > 5. When 
3.339 ^ 5^ < 5, there exists a 1-paramotor set of models with finite radii and 
a finite number, depending on n, with infinite radii. 
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1. INTRODUCTION 

This paper is the second in a series deahng with general relativistic star 
models. Here we consider static spherically symmetric models. The line 
element for these models can be written as 

ds'' = -e''^^^^dt^ + r{Xf\N{Xfd\^ + dn^] , (1) 

with 

rfO^ = de^ + sin^ 0dip^ , (2) 

where (?i(A) is the gravitational potential, r(A) is the usual Schwarzschild 
radial parameter, and A'^(A) a dimcnsionless (under scale-transformations) 
freely specifiable function. The choice N = 1 corresponds to isotropic 
coordinates and the function N can hence be viewed as a relative gauge 
function with respect to the isotropic gauge. The coordinate A is a spatial 
radial variable, defined by the choice of N. 

The matter content of the star is assumed to be a perfect fluid, described 
by the energy-momentum tensor 

Tab = PUaUb + P {gab + UaUb) , (3) 

where p is the energy density, p the pressure, and m° the 4-vclocity of the 
fluid. In this paper we focus on polytropic equations of state 

p = fc/ , (4) 

where the constant F is related to the polytropic index n according to 
r = l + l/n. It will be assumed that n > 0, where the case n = 
corresponds to an incompressible fluid. In the limit n ^ oo, the polytropic 
equation of state (4) becomes linear and scale-invariant (it thus corresponds 
to the special case po = in the linear equation of state, p = po + {r]— l)p, 
treated in [13]). 

Newtonian polytropic models have been studied for over a hundred years. 
Early results have been extensively described by Chandrasekhar [3], how- 
ever, it is still an active area of research as exemplified by the fairly recent 
papers by Kimura [9] and Horedt [8]. This indicates that the Newtonian 
case is quite non-trivial, and the relativistic case turns out to be even more 
complicated, as we will demonstrate in this paper. The history of the rela- 
tivistic case is also quite long. For example, Tooper [16] studied these mod- 
els 35 years ago, but they have also been studied recently by, for example, 
deFelicc et al [4] . Often polytropes describe low or high pressure regimes of 
more realistic equations of state for white dwarfs and neutron stars. These 
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regimes decide many of the models physical features, which therefore can 
be understood if the corresponding polytropic models are understood. One 
can also constnict physically realistic composite equations of state by let- 
ting certain pressure regimes be described by polytropic equations of state. 
There are thus ample reasons to study models with polytropic equations of 
state, and it is not surprising that these models have attracted attention 
for a long time. 

Substituting (1) into the equations of motion of the fluid, WaT"'^ = 0, 
relates the gravitational potential (p in (1) to the matter content of the 
spacetime according to 

# _ 1 

dp p + p 

To facilitate the study of the gravitational field equations, we introduce 
various sets of bounded dimensionless variables (under scale-transformations), 
which regularize the gravitational field equations. Recasting the equations 
into regularized form on a compact state space allows us to numerically 
and qualitatively explore the solution space, using methods from dynami- 
cal systems theory. It is worth noticing that one can modify the approaches 
we introduce and apply them to large classes of equations of state. The 
dynamical systems are 3-dimensional in all cases. This is a great advan- 
tage since one therefore can visualize the compact state space and obtain 
a picture of the structure of the entire solution space for large classes of 
equations of state. 

In the interior of a star, a mass function m(r) can be defined (see, for 
example, Misner & Sharp [11]). Buchdahl has derived inequalities limiting 
the behavior of the mass for solutions with a regular center 

m > ^nr^p , (6a) 
m < ^TTr^pc , (6b) 

— H (67rr^p - 1) + 4:nr'^pUnr'^p - 2) , (6c) 

r / r 

where pc is the energy density at the center of the star (see, for example, 
Buchdahl [1] and Hartle [7]). Alternatively, (6c) can be written as 

— < I (l - 67rr^p -I- a/1 + 67rr2p^ . (7) 

These inequalities are satisfied by all regular solutions, but are equalities 
only for the incompressible fluid. For such models, (6a) and (6b) charac- 
terize the regular solutions, while (6c) corresponds to a certain non-regular 
solution with a positive mass singularity. 



(5) 
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A static spherically symmetric perfect fluid model has to be isolated, 
i.e., it has to have a boundary at a finite radius, if it is to be considered as 
a star model. At the radius where the pressure of the fluid vanishes, the 
interior solution can be matched with the static Schwarzschild solution 

ds' = - (1 - ^) dt' + -^^^ + r'dn' , (8) 

(see Schwarzschild [15]). 
If (7) is evaluated at the surface of the star, defined by p = 0, we obtain 

¥ < § . (9) 

where M is the mass and R the radius of the star. 

The outline of the paper is as follows: Section 2 constitutes the main 
part of the paper. In this Section, the field equations are recast into two 
complementary regular dynamical systems on compact state spaces. Both 
systems are studied numerically and qualitatively. The results are com- 
pared and various physical implications are discussed. For example, it is 
shown that general relativistic models have finite radii when the polytropic 
index n satisfies < n < 3.339 and infinite radii when n > 5. When 
3.339 < n < 5, there exists a 1-parameter set of models with finite radii 
and a finite number, depending on n, with infinite radii. The situation 
thus differs from the Newtonian case where regular models have finite radii 
when the polytropic index satisfies < n < 5, otherwise not. In Section 
3 we discuss how one can adapt the different approaches of this paper to 
facilitate studies of more general barotropic equations of state. We con- 
clude in Section 4 with some remarks. In Appendix A, we indicate how 
the various variables were found. In Appendix B we show the relationship 
between our variables and the so-called Lane-Emden variables. 

Throughout the paper, geometric units with c = G = 1 are used, where 
c is the speed of light and G the gravitational constant. Roman indices, 
a, b, ... = 0, 1, 2, 3 denote spacetime indices. 

2. DYNAMICAL SYSTEMS AND DYNAMICAL SYSTEMS 

ANALYSIS 

We now present two complementary formulations adapted to the poly- 
tropic equation of state (4). In the first formulation, we introduce dimen- 
sionless variables 



{^,K,y} 



(10) 
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according to 

N'=y'K, ct>' = ^y, = ^ 2/=^' (") 

Stta \ y J P + P 

where the prime denotes differentiation with respect to the independent di- 
mensionless spatial variable defined by the above choice of N. The quantity 
P is defined by 

P = 1 - - , (12) 

and can be expressed as 

P = SirKr^p . (13) 

Thus, the assumption of a non-negative pressure implies P > 0. The 
gauge, however, breaks down when P = 0. Prom (11), it follows that K 
is positive. We assume that the energy density p and the pressure p arc 
both non-negative. This, together with the defining equation for y in (11), 
implies that the variable y satisfies < y < 1. 

For a non-negative polytropic index, n > 0, and for a positive constant 
fc in p = A:p^+^/", the subset y = corresponds to the limit p,p 0, while 
the subset y = 1 corresponds to the limit p,p — > oo. However, one can 
also vary k. Setting k = ^ p = 0, while k~^ = => p = 0. These 
limits corresponds to y = and y = 1, respectively, as follows from the 
definition of y in (11). The same "equations of state" (p = or p = 0) 
can be obtained from the linear equation of state, p = (7 — l)/o; p = 
is obtained when 7=1 and p = is obtained when 7 ^ 00. Thus the 
subsets y = and y = 1 can be viewed as state spaces associated with the 
linear equation of state, p = (7 — l)p, in the limits 7^1 and 7 — > 00, 
respectively. 

Integrating (5), leads to 

e'^ = a(l-j/)^+" , (14) 

where a is a freely specifiable constant corresponding to the freedom of 
scaling the time coordinate t in the line element (1). This, in turn, re- 
flects the freedom in specifying the value of the gravitational potential </> 
at some particular value of r. Matching an interior solution with the ex- 
terior Schwarzschild solution when p = 0, however, fixes this constant to 

a = y'l — It is worth noting the simple relationship between the 

variable S and the logarithmic derivative of the gravitational potential. 



# _ E 
d\nr " 1 - E 



(15) 
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In terms of the variables {'E,,K,u}, the gravitational field equations take 
the form 

E' = -yKT.+ ^P{l + 2y^AyY.) , (16a) 
K' = 2y{Y? - P)K , (16b) 
y' = -(l-r-i)y(l-j/)E , (16c) 

where T = 1 + i. It follows from (16b) and (16c) that X = and y = 0, 1 
are invariant subsets. By differentiating (13), we obtain 

P'= [-E(l + 22/ + 2zjE) + 2ify]P , (17) 

which shows that P — is an invariant subset as well. All of these invariant 
subsets constitute the boundary of the physical state space, and by includ- 
ing them in the analysis we obtain a compact state space. For the K — Q 
subset, the remaining equations describe plane-symmetric polytropic mod- 
els (compare with the linear equation of state discussion in [13]), and we 
therefore refer to this boundary subset as the plane-symmetric boundary. 
The state space is shown in figure 1, along with the different boundary 
subsets. 




Plane-symmetric boundary 
(K^O) 



FIG. 1. The state space for polytropic models, together with the different boundary 
subsets, in terms of the variables {'S,y,K}. 



Differentiating r in (11), yields 

r' = (1 - E)yr , 



(18) 
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which shows that r is a monotone function for y ^ 0, r ^ and S ^ 1. 
Together with monotone functions on the boundary subsets, this imphes 
that all attractors arc equilibrium points on the boundary subsets y = 0,1 
and P = 0. The equilibrium points of (16a)-(16c) and their corresponding 
eigenvalues are given in Table 1. See also figure 1. 

TABLE 1. 

The equilibrium points and their stability using the variable set {E, K,y}. 



Eq point 


E 


K 


V 


Eigenvalues 


Li 


Es 







-Ss , , 







1 


2/c 


, -yc , 2yc 


Pi 


2 


1 


1 


12 2 




3 


9 




3 ' 3 ' 3(l+n) 


P2 


1 





1 


1 2 

' ^ ' 1+n 


Ps 


-1 





1 


' ' ^ ' l+n 


Pa 


3 
4 





1 


7 1 3 

8 ' 4 ' 4(l+n) 



The solution structure of the boundary subsets is shown in figure 2. 

The equilibrium point Pi is associated with the self-similar Tolman solu- 
tion, discussed in [13], in the limit 7 — » oo when p = (7 — l)p. We therefore 
refer to this point as the Tolman point. A single orbit enters the interior 
state space from Pi (it is associated with the single positive eigenvalue of 
Pi). We refer to this orbit as the Tolman orbit, although the corresponding 
exact solution is not known (except in the incompressible fluid case). An 
approximation of this solution near the center has been given by deFelice 
et al [4] . There originates a 2-parameter set of orbits into the interior state 
space from the hyperbolic source P2. These orbits correspond to solutions 
with a negative mass singularity. There are no interior orbits associated 
with the hyperbolic saddle P3. A 1-parameter set of orbits enters the in- 
terior state space from the hyperbolic saddle P4. This point is associated 
with a self-similar plane symmetric solution discussed in [13] in the limit 
7 ^ c» when p = (7 — l)p. The orbits coming from this point are associated 
with a negative mass singularity. 

The Ss < part of the line Pi constitutes a transversally hyperbolic 
source from which a 2-parameter set of orbits originate. Solutions associ- 
ated with these orbits start out with a negative mass. Eventually, for a 
sufficiently large radius r, all solutions starting out with negative mass ac- 
quire a positive mass, since the pressure is increasing with r when the mass 
m(r) is negative and since the pressure and energy density are assumed to 
be positive. The Sg > part of the line Pi constitutes the global sink. We 
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FIG. 2. Orbits in the boundary subsets using the variables {2,7^,1/} for (a) the 
boundary subset y = 0, (b) the boundary subset y = 1, (c) the plane-symmetric bound- 
ary K = 0. (d) The P = subset projected onto the K = plane. 



note that the point Eg = on Li is non-hyperbolic, with aU three eigenval- 
ues equal to zero. This makes it hard to analytically, as well as numerically, 
to resolve the dynamics in the vicinity of this point. One possibility to gain 
further information is to consider center-manifold theory, sec Can [2] for 
an introduction to the subject. We will not pursue this possibility here. 
Instead wc will rely on an indirect discussion concerning the models radii 
and masses and on a second formulation, introduced below, which to a 
large extent resolves this problem. 

Each point in the equilibrium set L2 corresponds to the flat Minkowski 
geometry written on spherically symmetric form (compare with the treat- 
ment in [13] where this solution was represented by a single equilibrium 
point). A single orbit, associated with the eigenvalue 2yc, enters the inte- 
rior state space from each point of L2 when yc > 0. The 1-parameter set 
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of all such orbits form a separatrix surface describing the regular subset of 
solutions, since the corresponding solutions all have regular centers. The 
individual solutions arc parametrized by yc. which in turn is determined 
by the value of the central energy density, pc, and central pressure, pc, 
according to tin — ^1 — —ttt, — ■ Demanding a causal star model, i.e., a 
model where the velocity of sound of the fluid is less than or equal to the 
speed of light, yields the following inequality for the initial values yc at the 
center of the star: 

Vc < ^Ar- ■ (19) 
^ 1 + 2n ^ ' 

Where an orbit ends on Li is intimately connected with the solutions 
total mass and radius. The ratio between the mass function m(r) and the 
radius r, is given by (as follows from the definition of m given by Misner 
& Sharp [11]) 

It follows from the above equation that the static condition < ^ is 
automatically satisfied in the interior state space. Note also that it follows 
that S > for solutions with positive mass. Evaluating (20) on Li, where 
K = Ks = l-i:'^, yields 

- = ^. (21) 

A linear analysis shows that orbits ending at Eg > on Li correspond to 
solutions with finite radii and masses. Orbits ending at the non-hyperbolic 
point Eg = on Li describe solutions with infinite radii and finite masses, 
or masses that approach infinity slower than r. 

The Buchdahl inequalities limit the possible behavior of the regular sub- 
set. Expressed in the present variables the first inequality, (6a), takes the 
form 

2y(3S(l - S) - P) > P , (22) 
where P is given in (12). For regular orbits, characterized by yc, (6b) yields 

\ l+n y \ 1+n 

^) (l^) ■ '''' 

The third Buchdahl inequality, (6c), leads to 

K >{1- 2S)2 . (24) 
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On the surface this inequahty results in Sg < |. 

In figure 3a,b,c some orbits in the regular subset are shown, together 
with the important Tolman orbit. As seen from the figures, the regular 
orbit on the y = 1 boundary, connecting the y = 1 end of L2 and the 
Tolman equilibrium point Pi, and the Tolman orbit play important roles 
for understanding the regular solutions, the latter acting not only as a 
boundary but also as a "skeleton" orbit for the regular set, when n > 0. The 
situation strongly resembles that of the linear equation of state discussed in 
[13]. This is not a coincidence. The regular orbit on the y = 1 boimdary and 
the Tolman orbit are completely analogous to the corresponding solutions 
in the linear equation of state case. As discussed earlier the j/ = 1 set 
corresponds to the limit 7 ^ 00 when p = (7 — l)p. From the figures we 
also see that as n increases, the orbits are pushed towards the point Sg = 
on Li and that the maximum value of Eg for the regular solutions is closely 
connected with the value of Eg for the Tolman orbit. 

For the incompressible fluid case, n = 0, the Tolman orbit describes a 
"simple" lower boundary of the regular subset and gives an upper bound 
Eg = I for the regular solutions when they end at Li, see figure 3a. This 
upper bound is just the Buchdahl inequality (6c) , which holds for all equa- 
tions of state. When < n < 5, the situation is more complicated. In this 
case some regular orbits spiral around the Tolman orbit and some regular 
orbits end at larger Eg-values on Li than the Tolman orbit. Nevertheless, 
the largest Eg-value is close to the one determined by the Tolman orbit (for 
which Eg < I when n ^ 0), see figure 3b. When < n < 3 all orbits end 
on Eg > 0; when 3 < n < 5 a 1-parameter set of orbits end on Eg > and 
a finite number at Eg = 0, as follows from the results obtained from the 
second formulation below. When n > 5, all regular orbits spiral around 
the Tolman orbit and they all end at Eg = 0, see figure 3c. 

Combining the above discussion with equation (15) tells us that the 
behavior of the gravitational potential close to the surface of a regular 
solution (which might be located at infinity), is to a considerable degree 
limited by the behavior of the gravitational potential of the non-regular 
solution that corresponds to the Tolman orbit. Moreover, this behavior is 
intimately connected with the equation of state. 

An advantage of the above formulation is the regularity of the equations 
which makes the analysis of the center of the stars simple (compare with 
Rendall & Schmidt's analysis of singular differential equations in [14]). In 
addition the above formulation works well when the energy density is large, 
even though one might be far from the center. The formulation also works 
well for any values of the energy density when n is small. Unfortunately, it 
is not well suited to deal with the low energy density limit when n is close 
to 3 or higher. 
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FIG. 3. Orbits belonging to the regular subset using the variables {S, K, j/} for (a) 
n. = (b) < n < 5 (c) n > 5. 



To facilitate a study of polytropic stars at low energy densities for any 
non-negative value n, we now present another formulation of the polytropic 
equations. We again choose 



V = 



P + P 



(25) 



as a variable, but make a non-linear transformation of K and E to a new set 
of dimensionless variables U and V (see Appendix A). We first introduce 
the Newtonian homology invariants u,v (see, for example, Kimura [9]) 



Airr p mp 

u — , V — 

m rp 



and then define the variables U and V according to 



U 



A-Kr p 



1 + u 47rr^p + m/r 



m,/r 



1 + V m/r + p/p 



(26) 



(27) 



The assumptions of positive energy density, pressure, and mass lead to the 
inequalities < U < 1 and < < 1. It follows that 



k"UV fl-y 

4^(i-c/)(i-y) \ ^ 

ef = a(l-y)^+" . 



(1-1 



(28a) 
(28b) 



We now choose a new independent variable such that the function N 
takes the form 



7V2 ^ (1 - yf{\ - V){1 - UfFr~ 



F={l-V){l~y)^2yV , (29) 
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where has been given in equation (28a). We thereby incorporate the 
assumptions of positive energy density, pressure, and mass in the gauge, 
which becomes ill-dcfincd when these requirements arc not satisfied. 

In terms of the variables {U, V, y}, the gravitational field equations, which 
are again completely regular, takes the form 

U' = U{l-U)[{l~y){3-m)F-r-^G] , (30a) 

V' = V{1 - V) [(1 - y){2U - 1)F + (1 - r-i)G] , (30b) 

y' = -(1 - r-')y{l - y)G , (30c) 

where 



G = V[{l^U)il-y)+yU] , (31) 

and where F = 1 + 1/n. The prime denotes differentiation with respect 
to the independent variable, defined by (29). By including the invariant 
subsets U = 0,1,V = 0,1. y = 0. 1 of (30a)-(30c) we obtain a compact state 
space. These subsets can be interpreted in terms of limits of ^i. H = fcp^/" 

and pr^. The incorporation of the positive pressure, energy density and 
mass conditions in the gauge has led to that we now have invariant subsets 
corresponding to setting these quantities to zero. From (30c) it follows 
that y is a monotonically decreasing function, and all orbits in the interior 
of the state space approach the y = subset. The equilibrium points of 
(30a)- (30c) and their eigenvalues arc given in Table 2. 

Before continuing the discussion of the system (30a)- (30c), it is of interest 
to consider the Newtonian equations in terms of the variables {U,V,y}. 
These are given by: 

U' = J7(l-C/) [(3-4C/)(l-y)-r-V(l-Z7)] , (32a) 
V' = V{1-V)[{2U -1){1-V) + {1-T-'^)V{1-U)] , (32b) 
y' = _(i _ r-i)y(l - y)V{l - U) . (32c) 

Note that the non-homologous variable y is decoupled in the above system, 
leaving a coupled 2-dimensional reduced homology invariant set of equa- 
tions. The reduced homology-invariant set of equations is identical to the 
y = subset of the relativistic equations (30a) and (30b). We therefore 
refer to the y = subset as the Newtonian subset. Taking the Newtonian 
limit of all the relativistic equations (30a)- (30c) thus corresponds to letting 
y go to zciro (except in the equation for y where one must be more carehil) . 
Hence we obtain the Newtonian homology-invariant equations as the "low 
pressure" boundary of the relativistic state space. It is worth noting that 
the projection of the Newtonian orbits onto the subset y = 0, coincides with 
the orbits on the same subset, because of the decoupling of the y-equation. 
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We now return to the relativistic case. We are considering static models 
and thus the inequahty ^ < 1 must be satisfied. Solving for ™ in equation 
(27) leads to the expression 



m / V 



r \1-VJ \l-y 
Inserting this into ^ < 1 yields the inequality 

1-2/ 



(33) 



1 + 2/ 



> V . (34) 



The state space, together with the intersection of the "static surface" (1 + 
y)V — (1 — y) = with the {7 = and U = 1 boundary subsets, is depicted 
in figure 4. Note that the static surface "forces" the static solutions away 
from the line L5 and the point P3 , which thus are of little interest for the 
static solutions. 

Expressed in the present variables, the first Buchdahl inequality (6a) 
leads to 

U<1. (35) 

This inequality "forces" the regular subset away from P5. Inserting equa- 
tions (33) ami6TTr^p=^UVy^il-U)-'^il-V)-^il-y)-'^ into (6b) leads 
to a complicated expression, which we will not give. For regular orbits 
characterized by yc, equation (6b) yields 

A disadvantage with the present formulation is that the entire y = 1 
subset, in the {S, if, y}-formulation (see figure 2(b)), has been "crushed" 
into the line L3. We also note from Table 2 that this line is non- hyperbolic, 
with three zero eigenvalues. This means that studying the dynamics in a 
neighborhood of this line can be quite complicated. By relying on results 
from the {E, K, yj-formulation, wc can circumvent most of these problems. 
Although the previous formulation gives a better picture of what happens at 
largo energy densities, the present formulation gives a better picture of the 
low energy density regime. The two formulations arc thus complementary. 

When discussing the solution structure using the {S, /sT, y}-formulation, 
we emphasized the importance of certain special orbits, for example, the 
Tolman orbit. Thus it is important to identify this orbit in the present for- 
mulation. The first step in order to do this is to identify where the Tolman 
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^3 




FIG. 4. The structure of the state space using the variables {U,V,y}. The curves 
indicated on the U = and U = 1 subsets are the intersections between these subsets 
and the "static surface" (1 + y)V — (1 — y) = 0. 



TABLE 2. 

The equilibrium points and their stability using the variables {U,V,y}. 
Here A3 = - -, "7/, „ , and a = 1 + 22n - 7n^ 

(n — 2}(l + 3n) 



Eq. point 


U 


V 


y 


Eigenvalues 


Restrictions 


Li 


1 





yo 


(l~yo)' , , 




L2 


3 
4 





Vc 


-1(1 -yc)^ , \[l-ycf , 




L3 


Uo 





1 


0,0,0 




u 








yo 


3(1 -yo)' , -(1-yo)' , 




Ls 





Vb 


1 


0,0,0 




Le 


Uo 


1 





, -{1-Uo) , -(1-f/o) 


n = 


Pi 





1 + n 





n-3 1 1 






2+n 


2+n ' 2+n ' 2+n 




P2 





1 





nil 
1+n ' 1+n ' 1+n 




Pz 


1 


1 


1 


n 11 






1+n ' 1+n ' 1+n 




Pi 


n — 3 


2(l + n) 





(5 - n ± yS) , A3 


n > 3 


2{n-2) 


l + 3n 


P5 


1 


1 





0,0,0 





point from which it originates is located. This point corresponds to a point 
characterized by J7 = (3 + n)/2{2 + n) and V = m. the equihbrium set 
L3, as can be seen using the variable transformations given in Appendix 
A. From these transformations one also sees that the "plane-symmetric" 
point P4, in the {E, K, yj-formulation, corresponds to the equilibrium point 
U — Q,V — Q,y — 1, i.e., the intersection of the three lines L3, L4 and L5. 
We know from the previous analysis that this latter point is associated 
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with a 1-parameter set of solutions originating from it and that these solu- 
tions start out with negative mass. Transforming from the {S, i<'}-picture 
to the {U, y}-picturc shows that the 1-paramctcr set of orbits starts from 
U = 0,V = 0,y = 1 entering the "negative mass part" of the state space 
outside the cube. They eventually end at Li when the total mass becomes 
zero. This line is just an artifact of the gauge choice. Using the previous 
formulation one can continue the solutions through Li into the the interior 
"positive mass part" of the cube. The situation is analogous to that en- 
countered in the diagonal homothetic approach for self-similar spherically 
symmetric models when one matches the spatially self-similar part with 
the timelike self-similar part, as described in Goliath et al [5], [6]. The 
2-parameter set of orbits associated with P2 in the {S, K, y}-formulation 
starts from L3 and enters into the negative mass regime of the state space 
outside of the cube, but eventually pass through Li into the interior state 
space, in the {U, V, y}-formulation. The 2-parametcr set of orbits associ- 
ated with the Sg < part of Li (in the {S, K, y}-formulation) start from 
Li (in the {U, V, y}-formulation) entering the exterior state space and even- 
tually pass through Li into the interior state space. 

In the {[/, V, y}-formulation, all regular solutions start from the line L2, 
and are again parametrized by yc- For n = 3, the equilibrium point P4 
enters the state space through the point Pi, changing the stability of Pi 
from a source to a saddle. As n increases, P4 first changes from a node to 
a focus and then, at n = 5, to a center in the Newtonian subset. When 
n > 5 P4 is a local sink of the full state space. The point P4 corresponds 
to a special singular Newtonian solution. 

To be able to determine if an orbit corresponds to a model with finite 
mass or radius, it is useful to consider equation (28a) and the following 
equations 

r' 
m' 

where F ^ (1-V)(1- y) - 2yV. 

The radius of a model is definitely infinite if r' does not approach zero as 
A — > -|-oo. The corresponding statement holds for the mass, m, as well. The 
variable y is monotonically decreasing and all orbits end at the Newtonian 
boundary subset. This, in conjunction with the stability of the equilibrium 
points given in Table 2, implies the following: Only orbits that end at the 
hyperbolic sink P2 (or the equilibrium set Lg for n = 0) correspond to star 
models with finite radii. These models also have finite mass. If an orbit 



{l-y){l-U)Fr , 
(1 - Vni - y) ' 
Anil - U)(l - V)3 \ l-y 



S—n 



(37a) 
(37b) 

(37c) 



16 



U. S. NILSSON AND C. UGGLA 



ends at one of the other possible equihbrium points, Pi or P4, a model 
with infinite radius is obtained. If it ends at Pi it yields a model with 
finite mass and if it ends at P4 the corresponding model has infinite mass. 
This follows from a linear analysis of equations (28a) and (37a)-(37c). Note 
that the above statements also hold for the Newtonian models, described 
by the y = subset and the decoupled y-cquation (32c). 

We now consider the Newtonian subset y = 0, since this subset is of 
great importance also for relativistic stars. 

2.1. The Newtonian subset 

Newtonian polytropic stars have been studied extensively in the litera- 
ture. Chandrasekhar [3] primarily used the Lane-Emden approach, while 
the homology-invariant variables {u, v} were used by, for example, Kimura 
[9] and Horcdt [8]. Both Kimura as well as Horedt used dynamical sys- 
tems techniques in their investigations, but since they did not use bounded 
variables, their state spaces were not compact. 

Prom these previous analyses, however, it is known that regular Newto- 
nian solutions exist for all n > 0, and that when < n < 5 they give rise to 
finite star models while n > 5 lead to infinite regular models. These results 
are easily obtained using the variables {U,V}. The solution structure of 
the Newtonian subset for different values of n is shown in figure 5. The 
orbit corresponding to the regular solution (called the E-solution in the lit- 
erature) starts at y = on L2, and we refer to it as the regular Newtonian 
orbit. As indicated in the discussion above, the equilibrium point where 
the orbit ends, depends on the value of n. For the solvable incompressible 
fiuid case n = 0, the regular orbit is characterized by C/ = |, and end on 
the corresponding point on Lq. For < n < 5, it ends at P2. Thus, for 
these values of n, the radius and mass of the star model are finite. We also 
note that for < n < 3 there exists a monotone function. 



{1-U){l-Vf ' ^^^^ 



which satisfies 

Z 



= 2nC/(l - F) + (3 - n)V{l - U). (39) 



For n = 5, however, the regular Newtonian orbit ends at Pi, and thus 
corresponds to a model with infinite radius but with finite mass. Note 
that for n = 5, the point P4 is a center. This is not just a feature of 
the linear analysis. The case n = 5 is exactly solvable, see, for example, 
Chandrasekhar [3], so there exists an integral, which in our variables takes 



GENERAL RELATIVISTIC STARS 



17 



the form 

D{1 - Uf/^{1 - Vf/^ = eVUVH , (40) 

with 

H = 2{1- U){W - 3f + 4UV{1 -V)- 3(1 - U){1 - Vf , (41) 

where D is a constant. Setting D = yields the regular orbit. Negative 
values of D parametrizes the closed orbits around F4. These orbits corre- 
spond to non-regular models with infinite radii. The orbits that end at P2 
are parametrized by the positive values of D. For n > 5 all orbits end at 
P4, giving rise to models with infinite radii and masses. 




(d) (e) (f) 



FIG. 5. Orbits in the Newtonian subset j/ = in terms of the variables {U, V} for 
(a) n = 0, (b) < n < 3, (c) 3 < n < (11 + 8V2)/7, (d) (11 + 8V2)/7 < n < 5, (e) 
n = 5, and (f) n > 5. 



2.2. Regular relativistic stars 

The dynamics of the relativistic case is much more complicated than in 
the Newtonian case. In contrast to the Newtonian case, a projection of the 
orbits onto the plane y = docs not coincide with the orbits in the y = 
subset, since the equation for y is not decoupled. The regular subset is a 
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1-parameter set of orbits starting from L2, defining a separatrix surface in 
state space. In the interior of state space, the Tolman orbit constitutes the 
boundary of this separatrix surface. Since y is a monotone fimction, all 
orbits end at the Newtonian boundary. We will now discuss the behavior 
of the orbits in the regular subset for different values of n. 

2.3. The incompressible fluid models 

The behavior of a relativistic incompressible fluid (n = 0) is quite rem- 
iniscent of the Newtonian case. All orbits in the regular subset lie in the 
plane U = j, and thus end at C/ = | on Le. The regular subset projected 
on the Newtonian subset y = is shown in figure 6a. 

We note that it is possible to solve the equations exactly for the incom- 
pressible fluid, just as in the Newtonian case (see, for example, Tooper [16] , 
who uses the relativistic Lane-Emden equation). 

2.4. Models with < n < 3 

For models with < n < 3, Table 2 shows that P2 is a hyperbolic sink. 
The fact that y is a monotone function that pushes all orbits down to 
the Newtonian subset, in conjunction with the existence of the monotone 
function (38) in the Newtonian subset for < n < 3, implies that P2 is, in 
fact, the global sink and all orbits end at this point. Therefore the radii 
and masses of all regular stars are finite when < n < 3. This is consistent 
with a theorem given by Makino [10], which in the present context states 
that the radius of a regular model is finite when 1 < n < 3 (see Theorem 
1, p 60 in [10]). 

The regular subset projected onto the plane y = is shown in figure 6b, 
where the regular Newtonian orbit is shown as a dashed line. Also shown 
is the Tolman orbit, starting from L^. Orbits characterized by an initial 
value of j/c close to unity, start out close to L3, and then follow the Tolman 
orbit closely. Note that the separatrix surface formed by the regular subset 
of orbits folds over itself. 

2.5. Models with 3 < n < 5 

For 3 < n < 5, the point P2 is still the only sink, and the majority of 
orbits will end at P2. Hence most regular orbits are expected to end there 
too, but there are other possibilities for where the orbits of the regular 
subset can end. There exists a 1-parameter set of orbits that end at Pi, 
and thus gives rise to a separatrix surface in state space. The correspond- 
ing solutions have finite masses but infinite radii. The boundary of the 
separatrix surface associated with Pi is described by an orbit in the U = 
submanifold, an orbit in the Newtonian subset y = connecting the point 
Pi with P4, and a relativistic orbit ending at P4. The solution correspond- 
ing to this latter orbit have both infinite mass and radius. If the "regular" 
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(a) (b) 
FIG. 6. Orbits belonging to the regular subset and the Tolman orbit using the 
variables {U,V,y} for (a) the incompressible fluid, n = 0, and (b) < n < 3. The 
Newtonian orbit corresponds to the dashed line. 



separatrix surface intersects this separatrix surface or its boundary orbit 
in the interior state space, there will be an orbit, ending at Pi or at P4, 
respectively, instead of P2. The two surfaces could, of course, intersect each 
other several times, resulting in several infinite solutions. What actually 
happens depends on the value of n, as discussed below. 

We have been forced to rely on numerical simulations in order to decide 
where the orbits of the regular subset actually end. To systematically 
explore if there are regular orbits that end at Pi or P4, we proceed as 
follows. We first note that when 3 < n < 5, there always exist a constant 
Ci > 0, such that all orbits with yc < C\ end at Pi- This follows from 
our numerical investigation, but is also a consequence of Theorem 4 on 
p 994 in Rcndall & Schmidt [14] (sec also Theorem 2 p 64 in Makino 
[10]), which states the following in the present context. Any relativistic 
polytropic model with 1 < n < 5 and central density pc has finite radius 

if k is sufficiently small. Since yc = ^r+p~ ~ ^'^f^-j^x ' means that 
models with sufficiently small lead to finite stars. Since y is a monotone 
function, the intersections of the separatrix surfaces with a plane y = C2 
are curves. Hence, an intersection of the separatrix surfaces corresponds to 
an intersection of the two corresponding curves in the plane y = C2 > ci. 
To visualize the intersections of the two separatrix surfaces for a given n, 
we thus first ensure, numerically, that all orbits with yc < C2, where C2 > 
is chosen to be some convenient value, end at P2. We then numerically 
calculate the separatrix surfaces and plot their intersection with the plane 
y = C2. 

Numerical simulations indicate that all regular orbits end at P2 for 3 < 
n < 3.339, and all regular models are thus finite for this range of n, see 
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figure 7a. The first time the two separatrix surfaces intersect each other is 
when n « 3.339, see figure 7b. Thus there is a single regular orbit ending 
at Pi, corresponding to a solution with finite mass and infinite radius. All 
other regular orbits end at P2, giving rise to models with finite radii and 
masses. When n is larger than the above value there is always at least one 
separatrix intersection with a "Pi -orbit" and there is at least one solution 
with finite mass and infinite radius. 

In figure 7c the intersection between the two separatrix surfaces for n = 
3.6 with the plane y = C2 = 0.1 is shown. The single point of intersection of 
the two curves shows that there exists a single infinite solution. For n = 4.1, 
the result is shown in figure 7d. We note that there are several infinite 
solutions, but they arc all isolated, and this is typical. Intersections always 
seem to take place at isolated points. The behavior of the regular subset 
projected onto the Newtonian subset is shown in figure 8a for 3 < n < 3.339 
and in figure 8b for 3.339 ^ n < 5. One should note that as n approaches 
5, the behavior of the regular subset becomes more and more complicated. 

The boundary orbit of the Pi-separatrix surface, corresponding to the 
P4 orbit, may also intersect the regular separatrix surface. This happens 
when n 3.357 and n ~ 4.414. In these cases one obtains solutions with 
infinite masses and radii. 

The Tolman orbit, i.e., the boundary of the regular subset, also in- 
tersects the Pi-separatrix surface. This happens at an increasing rate 
as n approaches 5. The first five values of n when this happens are 
3.673, 3.939, 4.105, 4.221 and 4.309. 

The lowest initial value of yc leading to an infinite model decreases to- 
wards zero when n increases towards 5. For n = 5 it coincides with y^, = 0. 
The rclativistic 3-dimcnsional state space thus sheds light on the appear- 
ance of an infinite Newtonian star for n = 5. 

It is worth noting that the appearance of infinite solutions is a source of 
considerable problems in the Lane-Emden approach. 

2.6. Models with n = 5 

When n = 5, the separatrix surface that ends at the equilibrium point Pi 
completely encloses the regular subset of orbits. This makes it impossible 
for any regular orbit to end at Pj. The Newtonian regular orbit is the lower 
boundary of the separatrix surface and thus end at Pi , while all rclativistic 
regular orbits are forced towards P4 and the closed orbits surrounding P4, 
in the Newtonian subset y = 0. All rclativistic regular models therefore 
have infinite radii and masses. The enclosing separatrix surface is shown 
in figure 9a, and the projection of the regular subset and the Tolman orbit 
onto the Newtonian subset y = is shown in figure 9b. 
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u u 

(c) (d) 
FIG. 7. The intersection of regular subset and the separatrix surface ending at the 
equilibrium point Pi in the plane y = 0.1, for (a) n = 3.1, (b) n « 3.339, (c) n = 3.6, 
and (d) n = 4.1. 




2.7. Models with n > 5 

When n > 5, the situation is similar to the case n = 5. The separatrix 

surface associated with Pi now encloses all regular orbits (including the 
regular Newtonian orbit) and all regular orbits end at the sink P4. Hence 
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FIG. 9. (a) The behavior of the separatrix surface ending at the equilibrium point 
Pi that excludes any finite polytropic stars with n = 5. (b) Orbits in the regular subset 
and the Tolman orbit projected onto the Newtonian subset y = 0. 



all polytropic models with n > 5, Newtonian and relativistic, have infinite 
radii and masses. The enclosing separatrix surface associated with Pi is 
shown in figure 10a, and the projection of the regular subset and the Tolman 
orbit onto y = is shown in figure 10b. 




u 



W (b) 

FIG. 10. (a) The behavior of the separatrix surface ending at the equilibrium point 
Pi that excludes any finite polytropic stars with n > 5. (b) The orbits of the regular 
subset and the Tolman orbit projected onto the Newtonian subset y = 0. 



3. COMMENTS ON THE GENERAL BAROTROPIC 
EQUATION OF STATE 
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TABLE 3. 

The functional form of 1 — 1"^^ for some common equations of state. 
Here pm is the rest mass density and f) = T/{1 —T). 



Type Eq. of state 1 - T'^ 



Incompressible p = const 1 

Scale-invariant p oc p 

Polytrope p oc p^+i/" (1 + n)-^ 

Linear p={r]-l)p + po (1 - ??2/)/(l - ?/) 



Relativistic polytrope p = kp^ , p = pm+p/{T-l) {1 - r)y)/[n{l - y)] 



In this section we address the issue whether it is possible to use, or at 
least modify, the above formulations, in order to understand models with 
other barotropic equations of state. 

In general, a barotropic equation of state, p = p{p), leads to that 

r = ^ (42) 

is a function of y. For a given T{y), one can parametrize the equation of 
state in terms of y as follows: 



p = poexp(^J -j^:^^ — '^^^ —] , (43a) 



dy 

{T{y) - l)y{l - y) 



The parametrization, however, might break down, but in such cases one 
can in principle obtain the equation of state in a piece-wise manner, except 
when r = 1 everywhere. This latter situation corresponds to a linear scale- 
invariant equation of state, p = (7 — 1)/?, easily treated with other methods, 
see, for example [13]. 

To obtain a formulation for non-polytropic equations of state, one has 
to replace the constant F in the equations in the previous section with the 
function T{y). The functional form of the factor 1 — for some common 
equations of state is given in Table 3. 

It should be pointed out that in general y is not a monotone single- 
valued function of p. Several values of p and p might correspond to the 
same value of y. This can be visualized as follows: First note that there is 
a one-to-one correspondence between y and | since y = l_^_^p/p) ■ Then plot 
the equation of state in a p-p diagram. The slope of a straight line from 
the origin determines ^ and therefore y. From this construction one can 
deduce the functional properties of y as a function of p or p. If a straight 
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line intersects the equation of state curve several times, then several values 
of p and p correspond to the same value of y. This feature is of course 
encoded in the above equations. If this situation occurs, then y is not a 
monotone function and 1 — T(y)~^ becomes zero for one or several values 
of y in the interval < y < 1. For these values of y one obtains invariant 
subsets described by the scale-invariant linear equation of state. If there 
are only a few zeroes of 1 — F"^, one may match orbits and obtain the 
complete solution. 

It is only possible to obtain F as a function of y explicitly for a limited 
set of equations of state. Because of these difficulties, a better strategy 
for most equations of state, particularly if 1 — F"^ has several zeroes, is to 
change from y to another variable z = z{y), that is a monotone fimction of 
p, and treat y as a function of z in the equations and if necessary choose 
a new independent variable. Simple examples of new variables are, for 
example, p, some simple function of e"^ (like the redshift), or, if one wants 
a bounded variable, for example, z = -^rf^, z = J^j^^^ , where a? is some 
arbitrary constant with suitable dimension, but preferably associated with 
the equation of state of interest. It might happen that one obtains new 
boundaries replacing, for example, the y = 1 boundary for the physical 
interesting part. This is exemplified by the linear equation of state, for 
which 1 — F"'^ — ^i^y and the relativistic polytrope, for which 1 — F~^ = 
(1 — fjy)/fj{l — y). In these cases one obtains a boundary a,t y = ^ and 
2/ = 4 respectively, instead of at y = 1. Note that for the system to be 
useful, one requires that the equations are differentiable of at least order 
on the entire physical state space and its boundaries. 

To exemplify that one can use the polytropic formulations for other 
equations of state, we will first consider the linear equation of state and 
make comparisons with the approach described in [13]. In the {S,JC, j/}- 
formulation one only needs to replace 1 — y with 1 — rjy in the equation for 
y. Hence y = r]~^ < 1 replaces the boundary subset y = 1 and the phys- 
ically interesting region is characterized by y < 77"^ < 1. The resulting 
dynamical system is completely regular on the physical state space and its 
boundaries. 

In [13] the state space allowed orbits corresponding to solutions with 
negative pressure. This is not the case now. We have incorporated the 
positive pressure condition in the gauge and we therefore now have an 
invariant subset corresponding to setting p to zero. Instead of ending at a 
non- invariant surface of zero pressure, as in [13], the regular orbits now end 
at SgS > on Li, since all regular solutions have finite radii and masses, see 
[13]. In addition, all regular solutions now start from a line of equilibrium 
points L2, while in [13] they all started from a single isolated equilibrium 
point. 
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For a linear equation of state the {U, V, yj-formulation yields the follow- 
ing equations: 

U' = U{l-U)[{l-y){3-4U)F-r-^G] , (44a) 

V' = V{l-V)[{l-y){2U-l)F + {l^r-^)G] , (44b) 

y' = -y{l - rjy)G , (44c) 

where 

F = {l-V){l-y)-2yV , (45a) 

G = V[{l-U){l-y)+yU] , (45b) 

and where 

i_r-i = ^. (46) 

The expression {l — y)~^ occurs as a factor in the above dynamical system, 

except when 77 = 1 and = 0. This case corresponds to an incom- 
pressible fluid and can thus be obtained by setting n = in the polytropic 
equation of state. For r] > 1, y = r]~^ < 1 replaces y = 1 as the boundary 
subset of the physically interesting region. Thus y < ri~^ < 1, and this 
implies that I — y > 0. Hence the above system is completely regular on 
the physical state space and its boundaries. If one is so inclined, one can 
change the independent variable with a factor of (1 — y) so that one obtains 
polynomial equations. Note that y is a monotone function in the physi- 
cally interesting region for this case and thus there is no need to change to 
another variable z = z{y). 

Solutions that start with negative mass enter the interior state space 
through a line of equilibrium points, which did not exist in the linear equa- 
tion of state treatment in [13]. In addition, all regular solutions start from 
a line of equilibrium points, while they all started from a single isolated 
equilibrium point in [13]. Moreover, instead of ending at a surface of van- 
ishing pressure as in [13], they all end at the equilibrium point P2 on the 
y = subset. 

The pictures are thus quite different, but one can extract all interesting 
physical information about, for example, the regular solutions from any of 
the above formulations. However, note that it is easier to extract informa- 
tion about total masses and radii using the formulation given in [13]. 

The relativistic polytropic equation of state 

P = fcPm > P = Pm + ' (47) 
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where pm is the rest mass density and p the energy density, provides another 
example where the previous formulations works well. For this equation of 
state the polytropic exponent f coincides with the adiabatic index 

■■-(^)|. (-) 

which is not the case for F in p = kp^. This equation of state asymptot- 
ically approach the linear equation of state p — (T — l)p when p, p — !■ oo 
and the polytropic equation of state p = kfF when p, p 0. Hence, 
the y = 1 boundary is replaced with a boundary y = ^, described by a 
scale-invariant linear equation of state. The y = boundary is just the 
Newtonian polytropic boundary discussed previously. 



4. CONCLUDING REMARKS 

We have obtained two complementary formulations with regularized equa- 
tions on compact state spaces for polytropic equations of state. This has 
allowed us to obtain a global picture of the solution space of these models. 
The two approaches revealed that certain solutions play an important role 
for the remaining solution structure. 

The first approach (the {S, iiT, y}-formulation) showed that the regu- 
lar scale-invariant solution, the self-similar Tolman solution, and a special 
non-regular solution (corresponding to the Tolman orbit) to a large extent 
determine the high energy density behavior of the regular solutions. Note 
that the non-regular Tolman solution is scale-invariant towards the origin. 
An approximate expressions for the metric coefficients for this solution has 
been given by deFelice et al [4]. 

The second approach (the {U, V, y}-formulation) revealed that when n > 
3 there exists another special solution (the one ending at P4) of great 
importance for the remaining solution structure. When 3 < n < 5 other 
solutions spiral around this "skeleton" solution. When n > 5, the regular 
solutions approach this solution and end at its endpoint P4. This point 
is associated with a non-regular Newtonian solution, which describes the 
behavior of the regular solutions close to the surface when n > 5. 

The dynamical behavior is particularly complicated when 3 < n < 5. 
It is interesting to compare with Theorem 4 p 994 in Rendall & Schmidt 
[14], which in the present context states that regular relativistic polytropic 
models with 1 < n < 5 have finite radii if j/c is sufficiently small, and 
Theorem 1 p 60 in Makino [10], which in the present context states that 
when 1 < n < 3, the radius of a regular model is finite. Makino also 
comments that the restriction of Rendall & Sc;hnii(it may be incndtable when 
3 < n < 5 (p 62 in Makino [10]). For polytropes we have numerically shown 
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that all regular models are finite when n < 3.339. Thus, in this polytropic 
index range the restriction is not necessary. However, our investigation 
shows that for larger values of n it is. Our work shows that what is actually 
happening when 3 < n < 5 is determined in the relativistic regime and 
not the Newtonian regime, although this latter regime determines a finite 
number of possibilities. Considering the complicated dynamical behavior, 
it is hard to see how one could obtain our results except by numerics. 
Perhaps one could gain further insights by going beyond the static and 
spherical assumptions. For example, the behavior might be partly related 
to the stability features of the models. 

We have shown that the present methods can be generalized to other 
equations of state. The state spaces in such cases arc also 3-dimcnsional 
and this will make it possible to visualize the entire solution space for 
such equations of state as well. The first {S, /C, yj-formulation should be 
useful when an equation of state behaves linearly or polytropically when 
the energy density goes to infinity. For such equations of state one could 
expect that the regular scale-invariant solution, the self-similar Tolman so- 
lution, and a special non-regular solution (the Tolman orbit), determine, 
to a large extent, the high energy density behavior of the regular solutions. 
The second {J7, F, y}-formulation should be useful for equations of state 
that behaves polytropically when the energy density goes to zero. For such 
cases the low energy density regime should be described by the Newtonian 
boundary. An interesting question is what happens with the radius of a 
regular model when one has an equation of state that approaches a poly- 
trope with index 3 < n < 5 asymptotically when p, p ^ 0. It is probably 
determined in the relativistic regime, which depends on the actual equation 
of state. The relativistic polytropic equation of state nicely illustrates the 
above discussion. Since this equation of state is of considerable physical 
interest, we will treat this case in more detail in a future paper. 

We think that the present type of approach nicely complements those 
used by Rendall & Schmidt [14] and Makino [10] and that they together 
should constitute useful tools for further explorations of relativistic star 
models. 

APPENDIX A 

In this Appendix we will show how the different sets of variables can be 
derived. Starting with the line element in the form 



(A.l) 
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and introducing variables according to 

e = i) , a = 9 , B = e'^-'^ , (A.2) 

where a dot denotes differentiation with respect to t, leads to the following 
expressions for the gravitational field equations 

e = -26'^ + ea + B'^ + m-Kp , (A.3a) 

& = -2ea + cr^ + 47r(p + 3p) , (A.3b) 

B = {-9 + a)B , (A.3c) 

Snp = 9^-a^-B^ . (A.Sd) 

This system looks very much like systems in spatially homogeneous cosmol- 
ogy (see, for example, Wainwright & Ellis [17]), even though the physical 
interpretation is quite different. Hence one can import ideas from treat- 
ments of spatially homogeneous models to the present situation. 

To obtain the {S, K, j/}-formulation, we use similar variables to the so- 
called expansion-normalized variables used frequently in spatially homoge- 
neous cosmology (see Wainwright & Ellis [17]), and a dimensionless matter 
variable y according to 

S = ^ K = ^, y=J^, (A.4) 
9 9-^ p + p 9-^ 



The new independent variable is determined by 

dX 9 ■ 



(A.5) 



To obtain the {U, V, y}-formulation 2, we start with the Tolman-Oppenheimer- 
Volkoff approach 

dp ^ {p + p){m + 4TTr^p) (A 6a) 

dr r{r — 2m) ' 

dm ,9 / . „, N 

__ = 47rrV , (A.6b) 
dr 

(see, for example, [12]). We then introduce the variables {U,V,y} in 
terms of r,m,p, and p, according to the definitions (27) and (11), and a 
new independent variable defined by 

dv 

— = (l-y)(l-C/)F, F={l-V){l-y)-2yV . (A.7) 
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The relationship between the variables {U, V} and {S, K} is given by 



A.l. LANE-EMDEN VARIABLES 

In the Lane-Emden approach, dimensionless variables variables 0le, vle, 
and ^LE are introduced according to 



'LE — 



1/n 



A^m(r) ^ ^ 
^'LE = — ; , ?LE = Ar , 



(A.9) 



where ^4 is a constant defined by 



A: 



' 47r(l - yc)pc 
(1 + n)yc 



(AlO) 



(see, for example, Tooper [16]). 

The relationship to the variables {S, K, y} is given by 



7LE 



vle = 



;le 



yj^ - Vc) 
2/c(i - y) ' 



8(1 + n)3 
1 - - 

2if(i + n) Vi-y. 



1-y 



1+n 



1/2 



n-1 / 1 \ l+n1 -'-/^ 

1-2/ 



(A. 11a) 
W , (A.llb) 

, (A.llc) 



where 



(A12) 
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The relationship to the variables {U, V, y} is given by 
2/(1 - Vc) 



7LE 



;le 



yc(i-2/) ' 



(i + n)3(i-;7)(i-i/)3 Vyc(i-y) 



3 — n 



1/2 



uv 



Veil - y) 
{l + n){l-U){l-V) \yil-ye) 



n-ll 1/2 



(A. 13a) 



, (A. 13b) 



(A.13c) 
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